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abstract 


Both  exhaustive  and  non-linear  spectral  fitting  procedures  are 
applied  to  far -field  surface  wave  data  for  determining  seismic  source 
parameters,  especially  depth.  These  procedures  were  applied  to  a 
series  of  events  with  decreasing  control  on  the.  source  parameters. 

High  quality  Rayleigh  wave  data  from  the  Southeastern  Missouri  Earth¬ 
quake  of  October  21,  1965  was  analyzed,  and  the  parameter  corrobora¬ 
tion  with  independent  estimates  was  good,  with  excellent  spectral  fits. 
Two  low  magnitude  (mb  =  4.  1)  Italian  events  were  analyzed  individually 
and  by  taking  spectral  ratios.  Both  techniques  agreed  with  a  nearby 
bodywave  solution  in  the  same  tectonic  region.  Next,  six  events  from 
Sinkiang  were  analyzed,  again  using  both  techniques.  Results  were 
mixed,  but  with  better  azimuthal  coverage  and  more  realistic  effective 
Q  corrections,  spectral  fits  and  depth  estimates  of  high  quality  can  be 
expected.  Finally,  initial  efforts  are  directed  towards  determining 
the  importance  of  the  parameter  variations;  it  was  found  that  depth 
has  the  greatest  effect  on  the  shape  of  the  spectra. 
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ABSTRACT 


Both  exhaustive  and  non-linear  spectral-fitting  procedures  are 
applied  to  far-field  surface  wave  data  for  determining  seismic  source  para¬ 
meters,  especially  depth.  The  spectral  fit  solution  is  based  on  a  minimum 
residual  criterion  between  the  theoretical  and  observed  spectra.  The  initial 
efforts  are  directed  towards  determining  the  importance  of  the  parameter 
variations;  it  was  found  that  depth  has  the  greatest  effect  on  the  shape  of  the 
spectra.  We  also  attempted  to  fit  a  double  couple  source  to  theoretical  ex¬ 
plosion  data  using  only  Rayleigh  waves.  Although  our  analysis  yielded  an 
event  at  15  km  depth,  indicating  that  an  improper  fit  is  possible,  considera¬ 
tion  of  radiation  patterns  and  spectral  levels  demonstrate  that  this  is  an  ex¬ 
ceptional  case. 

The  automatic  spectral  fitting  procedures  were  applied  to  a 

series  of  events  with  decreasing  control  on  the  source  parameters.  This  pro- 

* 

cedure  was  supplemented  by  the  use  of's.  series  of  narrow  band  filters  for  re- 

\ 

moving  multipath  interference,  which  was  Successful.  High  quality  Rayleigh 
wave  data  from  the  Southeastern  Missouri  Earthquake  of  October  21,  1965 
was  analyzed,  and  the  parameter  corroboration  with  independent  estimates 
was  good,  with  excellent  spectral  fits.  Two  low  magnitude  (m,  =  4.  1)  Italian 
events  were  analyzed  individually  and  by  taking  spectral  ratios.  Both  tech¬ 
niques  agreed  with  a  nearby  bodywave  solution  in  the  same  tectonic  region. 
Finally,  six  events  from  Sinkiang  were  analyzed,  again  using  both  techniques. 
Results  were  mixed,  but  with  better  azimuthal  coverage  and  more  realistic 
effective  Q  corrections,  spectral  fits  and  depth  estimates  of  high  quality  can 
be  expected. 
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SECTION  I 


INTRODUCTION 

The  problem  of  obtaining  seismic  source  parameters  from 
surface  wave  radiation  patterns  has  been  investigated  by  several  authors  over 
thr-  past  few  years.  Lambert  et  al.  (1972)  combined  experience  in  using  surface 
wave  data  with  a  trial-and-error  approach  to  obtain  their  solutions.  Mitchell 
(1973)  used  partial  derivatives  of  the  source  parameters  to  guide  a  trial-and- 
error  approach  in  fitting  surface  wave  radiation  patterns  from  the  Southeastern 
Missouri  Earthquake  of  1965.  Tsai  and  Aki  (1970a,  b),  Tsai  (1972a,  b),  Tsai 
and  Shen  (1972),  and  Weidner  (1972)  used  essentially  an  exhaustive  comparison 
scheme  between  the  observed  and  theoretical  spectra,  with  the  model  which 
yielded  the  minimum  residual  chosen  as  the  solution.  Turnbull  and  Alexander 
(1973)  refined  this  approach,  using  a  non-linear  iterative  scheme  on  the 
parameters  dip( 5 ),  slip(A),  strike(</>),  and  seismic  moment(M)  combined  with 
an  exhaustive  variation  with  depth.  The  following  report  represents  an 
attempt  to  refine  the  work  of  Tsai  and  Turnbull  into  a  unified  approach,  and 
combine  the  strongest  parts  of  their  methodology. 

In  Section  II,  each  part  that  contributes  to  the  seismic  record 
is  discussed  in  detail.  Our  analysis  of  the  amplitude  snectra  of  surface 
waves  examines  the  energy  in  the  frequency  band  from  0.  02  to  0.  10  Hz  (T 
from  10  to  50  seconds).  Because  the  theoretical  spectra  of  the  Rayleigh  and 
Love  waves  are  calculated  from  a  layered  half  space  of  nondissipating  elastic 
medium,  corrections  must  be  made  for  attenuation  due  to  anelasticity.  A 
discussion  of  this  correction  is  given,  and  the  attenuation  for  a  particular 
source- station  path  is  given  in  Section  VI  in  the  analysis  of  the  Sinkiang 
events.  The  observed  spectra  also  are  affected  by  interference  due  to 


multipath  arrivals.  This  interference  takes  the  form  of  spectral  modulations. 
Their  removal  was  accomplished  using  a  series  of  narrow  band  fitters  as 
discussed  by  Alexander  (1963).  Finally,  i„  Section  II,  the  source  excitation 
function  is  discussed  in  terms  of  the  source  time  function,  the  amplitude 
response  of  the  medium,  the  radiation  pattern,  and  the  application  of  Tsai's  and 
Turnbull's  combined  methodology  to  obtain  seismic  source  parameters. 

This  report  is  essentially  a  discussion  of  analysis  techniques 
applied  to  surface  wave  data  with  various  degrees  of  control.  Section  III, 
though ,  represents  the  beginning  of  a  detailed  study  of  the  intrinsic  capabilities 
of  the  methodology.  If  we  do  not  have  independent  determinations  of  the  source 
parameters,  such  as  the  dip(5 ),  slip(A),  and  strike!*)  from  body  wave  data, 
and  depth  from  pP,  we  must  ask  ourselves  the  meaning  of  the  minimum 
residual  solution.  We  are  investigating  this  problem  on  two  fronts.  First, 
amplitude  spectra  and  radiation  patterns  were  generated  for  many  repre-  ’ 
sentative  cases.  From  these,  conclusions  can  be  reached  on  parameter 
variations  and  spectral  shapes.  Secondly,  in  applying  our  theoretical  double 
couple  source  model  to  the  data,  we  have  always  assumed  that  the  event  was 
an  earthquake.  What  if  the  data  were  really  generated  by  an  explosion?  A 
discussion  of  the  application  of  an  earthquake  source  to  explosion  data  is 
given  for  several  cases. 

The  next  three  sections  present  applications  of  the  methodology 
to  multisite  surface  wave  data  with  varying  degrees  of  control  on  the  source 
parameters.  Each  successive  section  analyzes  data  with  less  control  than 
the  previous  one.  In  this  way.  we  can  get  some  idea  of  the  reliability  of  our 
method.  In  Section  IV,  superior  quality  Rayleigh  wave  amplitude  spectra 
from  the  Southeastern  Missouri  Earthquake  of  1965  are  analyzed.  Prev.ous 
solutions  for  the  source  parameters  have  been  obtained  using  both  the  surface 
wave  and  body  wave  data,  so  comparison  with  our  result  is  possible.  This 
case  demonstrates  what  is  possible  when  extensive  multi-azimuth,  non- 
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multipathed  data  with  good  attenuation  correctiona  are  available.  We  also 
examine  the  effect  of  source  structure  on  our  solution. 

Section  V  presents  the  analysis  of  data  for  two  low  magnitude 
Italian  events  (mb  about  4.0)  where  only  three  stations  (VLPE)  were  recording. 
Each  event  is  analyzed  separately,  and  then  spectral  ratios  of  the  events  are 
used  to  obtain  solutions.  From  body  wave  solutions  for  events  in  the  same 
tectonic  region,  we  can  see  that  our  solution  is  reasonable. 

In  Section  VI,  the  analysis  of  six  events  from  the  Sinkiang 
region  of  southern  Asia  is  presented.  The  problems  inherent  in  de- 
multipathing  and  in  the  application  of  an  attenuation  correction  for  events 
from  this  region  recorded  at  VLPE  stations  are  discussed.  In  particular, 
an  intensive  effort  was  made  to  determine  the  attenuation  coefficient  for  the 
Sinkiang-Chiang  Mai  (CHG)  travel  path  using  explosion  data.  No  corroborative 
evidence  of  source  parameter  values  were  available  at  this  time. 

Finally,  in  Section  VII,  we  give  a  summary  of  the  results 
given  in  the  previous  sections,  and  the  paths  this  investigation  expects  to 
take  in  the  future. 
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SECTION  II 


A  DISCUSSION  OF  THE  TRANSFORMED  SEISMIC  RECORD  AND 
THE  APPLICATION  OF  SPECTRAL  FITTING  PROCEDURES  TO  THEM 

Using  an  exhaustive  approach  or  a  non-linear  regression 
analysis  on  synthetically  generated  data  as  a  test  for  the  method  can  be  quite 
deceptive.  It  is  necessary  as  a  test  for  the  correctness  of  the  software,  but 
virtually  useless  as  a  measure  of  its  practical  applicability.  Considering  the 
average  quality  of  the  VLPE  data  used  in  Section  V  and  VI  of  this  report, 
parameter  estimates  obtained  from  Rayleigh  and  Love  wave  data  recorded 
at  two  sites  are,  in  the  opinion  of  this  author,  at  best  a  first  guess.  The 
reasons  for  this  statement  lie  in  a  detailed  analysis  of  those  factors  which 
perturb  the  observed  data  to  such  a  degree  that  our  theoretical  description 
of  the  seismic  source  is  inadequate. 

We  can  begin  such  an  analysis  by  examining  each  of  the 
component  parts  of  the  Fourier  transform  of  the  seismic  record.  The  trans¬ 
form  of  the  record  obtained  at  station  i  resulting  from  event  £  can  be  written: 

R.  .  (o>)  =  I.(oj)  \3m  {oj,  r) 
if  i  l 

where:  I. (a*)  -  Instrument  response 

i 

U^  (  (d  ,  r)  -  Frequency  domain  displacement 

The  instrument  response  at  each  of  the  VLPE  sites  is 
furnished  on  regular  basis,  so  little  error  should  result  from  this  correction. 
The  frequency  domain  displacement  can  be  written  as 

U^  (co,r)  =  C £  (w,  r)  Sj^(oj)  e  ^ 

II-l 


where:  G^(co,r)  -  Path  transfer  function 

S^(co)  -  Source  excitation  function 
-  Combined  phase  terms 

The  path  transfer  function  is  probably  the  weakest  link  in  our  analysis.  It 
involves  essentially  two  corrections  to  the  observed  data;  a  correction  for 
effective  attenuation  and  elimination  of  multipathing.  Two  approaches  for 
obtaining  an  effective  Q  are  given  in  Appendix  A.  One  involves  two  or  more 
stations  nearly  on  the  same  great  circle  path  (amplitude  equalization),  and 
the  other  involves  multisite  data  from  an  explosion  yielding  an  average  value 
for  a  region  (Tryggvas son's  method).  In  both  methods,  high  quality,  multi¬ 
site  data  is  required.  At  present,  these  numbers  are  not  available  for 
critical  source- station  paths  from  the  Eurasian  continent,  and  a  concerted 
effort  to  obtain  these  values  is  highly  desirable. 

Elimination  of  multipathing  has  been  accomplished  by  the  use 
of  a  series  of  narrow  band  filters.  This  technique,  first  used  by  Alexander 
(1963)  on  surface  wave  data,  is  described  in  detail  in  Appendix  R.  The 
application  of  these  filters  produces  a  series  of  wave  packets  corresponding 
to  the  center  frequencies  of  each  filter.  (These  packets  are  used  in  the 
analysis  of  the  Italian  (Section  V]  and  Sinkiang  [Section  VI]  events).  If  multi¬ 
pathing  exists,  more  than  one  packet  will  occur  for  the  selected  frequencies. 

By  determining  the  group  velocity  of  the  maximum  amplitude  peak  of  each 
packet,  and  comparing  these  velocities  to  standard  group  velocity  curves, 
the  correct  packet  can  be  chosen.  The  maximum  amplitude  peak  of  this  packet 
is  then  the  value  of  the  spectral  amplitude  corresponding  to  the  packet 
frequency. 

There  are  three  main  sources  of  inaccuracies  in  this  type  of 
analysis.  If  leakage  occurs  between  the  components,  either  through  recording 
defects  or  improper  analytical  rotation,  the  results  are  essentially  meaningless. 
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A  general  shift  in  a„  the  packets  can  occur  front  group  delays  of  the  instru¬ 
ments  and  imprecise  location  of  the  event.  Careful  analysis,  though,  can 
minimize  all  of  the  above. 

This  technique  was  chosen  over  the  complex  cepstrum 
analysts  previously  applied  by  Tsai  ,1972a,  because  of  difficulties  in  inter¬ 
pretation.  The  complex  cepstrum  technique  (Linville,  1,71,  assumes  that  the 

r‘  7' *  Si8nal  h“  eSSentUUy  ‘he  —  «■»»«*■»  as  the  original  signal, 

n  act,  each  frequency  is  affected  differently  by  lateral  velocity  changes 

and  the  interpretation  of  the  cepstrum  becomes  highly  subjective.  The  ’ 

cepstrum  technique  seems  more  suited  for  multiple  event  analysis. 

Finally,  w.  can  write  the  source  excitation  function  as: 

S(  to)  =  M/  L(  to)  kt  to,  to,  h)  Xto, 

where.  -  Seismic  moment 

"  Rayleigh  or  Love  wave  number 
L/  (at)  -  Transform  of  source-time  function 
N/  M  -  Medium  amplitude  response 
X£{o))  -  Complex  radiation  pattern 

The  theoretical  earthquake  source  is  represented  by  a  double 
couple  which  varies  as  a  step  function  in  time  [source-time  function  La  <t  ,1 

whose  strength  is  M/,  the  seismic  moment.  A  detailed  discussion  cihe 
source  time  function  is  given  in  Appendix  C. 

The  medium  amplitude  response  N/  to)  represents  the  response 
of  the  source  structure  to  a  Dirac  delta  function  point  source.  This  structure 

”  m  ‘he  idMliZed  f°rm  °f “  Parallel  layered  half  space.  In  Anpendix 

’  SCVera‘  rePreSenli“i™  structures  are  displayed  with  their  corresponding 
Rayleigh  and  Love  wave  group  velocity  curves.  The  effect  of  structure  on 
the  determination  of  source  parameters  is  discussed  in  Section  IV  in  the 
analysis  of  the  Southeastern  Missouri  Earthquake  of  1,65. 
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The  determination  of  the  source  parameters  is  achieved  using 
a  combination  of  Tsai's  Method  for  a  localization  of  a  best  fit,  and  then  a 
non-linear  regression  scheme  using  the  localised  parameter  values  as 
initial  values  to  obtain  the  absolute  minimum  residual  solution.  A  detailed 
review  of  Tsai's  Method  is  given  in  Appendix  E,  while  the  mathematics 
involved  in  the  non-linear  regression  are  given  in  Appendix  F.  The  advantage 
of  using  this  form  of  analysis  is  that  it  reduces  by  a  significant  factor  the 
number  of  computer  runs  necessary  to  obtain  a  minimum  residua’  solution. 
Tsai's  Method  has  also  been  modified  in  terms  of  displaying  the  minimum 
residual  solutions.  Instead  of  examining  only  the  28  minimum  residuals 
which  correspond  to  28  combinations  of  dip  and  slip  angles,  distributions  of  the 
minimum  residual  for  each  parameter  over  all  the  parameter  combinations 
are  displayed.  Examples  of  these  displays  are  given  in  the  next  four  sections. 

The  need  of  an  expanded  use  of  the  minimum  residuals  lies 
in  the  quality  of  the  observed  data.  If  the  data  quality  is  poor,  knowledge  of 
the  distribution  of  the  minimum  residuals  can  help  in  making  the  decision  of 
whether  the  absolute  minimum  residual  is  really  the  best  parameter  solution. 

In  other  words,  the  distribution  of  the  minimum  residual  solutions  should  be 
used  in  making  a  choice  between  minimum  residuals  which  have  values  close 
to  each  other.  Also,  if  the  use  of  the  absolute  minimum  residual  value  of 
a  parameter  is  difficult  because  of  relatively  flat  distribution,  then  a  corres¬ 
ponding  degree  of  confidence  should  be  expressed  in  that  choice. 
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SECTION  III 


CORRELATION  BETWEEN  MINIMUM  RESIDUAL  SOLUTIONS 
AND  SPECTRAL  CHARACTERISTICS 

Before  we  examine  minimum  residual  solutions  from  data 
with  varying  degrees  of  control,  we  must  first  ask  ourselves  the  meaning  of 
the  minimum  residual  solution.  This  section  represents  the  beginning  of  a 
detailed  study  of  the  uniqueness  or  non-uniqueness  of  source  parameters 
(particularly  depth)  determined  from  far-field  surface  wave  data.  We  are 
approaching  this  analysis  asking  the  following  questions: 

1)  Since  the  minimum  residual  solution  represents  the  closest 
match  of  spectral  shape  between  the  theoretical  and  observed 
spectra,  what  source  parameter  variation  produces  the  greatest 
change  in  spectral  shape7 

2)  Is  there  any  readily  observable  spectral  property  which 
characterizes  a  source  parameter's  value,  such  as  a  spectral 
indicator  for  depth? 

3)  What  property  of  the  source  description  or  of  the  layered  half 
space  that  it  is  placed  in  affects  spectral  shape  the  most? 

4)  In  applying  a  theoretical  double  couple  source  description  to 
the  observed  data,  we  assume  that  the  event  is  an  earthquake. 
What  if  the  data  were  really  generated  by  an  explosion? 

In  this  section,  we  hope  to  supply  answers  to  questions  l,  2,  and  4,  with 
question  3  the  subject  of  future  studies. 

To  address  questions  l  and  2  spectra  and  radiation 
patterns  were  generated  using  a  theoretical  double  couple  point  source 
placed  in  a  Outenburg  plane  layered  half  space  (see  Appendix  F  for  the 
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theoretical  formulation).  Although  we  generated  spectra  and  radiation 
patterns  for  many  source  parameter  combinations,  a  representative  set  is 
displayed  in  Appendix  G.  From  this  set  of  figures,  we  can  make  the 
following  observations: 

•  Depth  variations  can  change  spectral  amplitudes  of  Rayleigh 

waves  over  three  orders  of  magnitude,  making  this  parameter, 
by  far,  the  most  sensitive  to  spectral  variation  (Figures  G-la, 

G ..5a,  C  -6a).  Love  wave  amplitude  spectra  have  a  smoother 
variation  over  two  orders  of  magnitude  (Figures  G-lb,  G-5b, 
G-6b).  Dip  variations  can  change  spectral  amplitudes  of 
Rayleigh  and  Love  waves  about  l.  5  and  l.  0  orders  of  magnitude 
respectively  (Figures  C,-2a,b,  G-7a,b).  Slip  variations  are 
about  the  same  as  that  of  the  dip  angle  (Figures  G-3a,  b, 

G-8a,  b),  with  strike  variations,  in  some  cases,  very  minimal 

(Figures  G-4a,  b,  G-9a,  b). 

•  For  a  vertical  strike- slip  fault,  Rayleigh  wave  amplitude 
spectra  exhibits  a  hole  which  changes  frequency  with  depth 
(Figure  G- 7a).  This  hole,  though,  is  greatly  reduced  for  a 
less  than  perfect  strike-slip  fault  (Figure  G-5a),  and  completely 
vanishes  for  a  dip- slip  fault  (Figure  G-6a). 

•  Rayleigh  spectral  variations  of  slip  for  a  vertical  fault 
exhibits  a  common  spectral  point  (Figure  G-3a)  which  moves 
across  the  spectrum  as  depth  varies  (Figures  not  shown). 

This  point  vanishes,  though,  for  non- vertical  faults  (Figure 

G-8a). 

These  observations  from  the  amplitude  spectra  indicate  that 
depth  is  the  most  sensitive  parameter,  but  no  one  spectral  indicator  for 
depth  is  readily  obvious. 
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Our  analysis  of  question  4  yielded  quite  interesting  results. 
Theoretical  spectra  for  a  100  kiloton  explosion  in  granite  at  shallow  depth 
(3  km)  was  generated  for  Rayleigh  waves  using  an  approximation  of  the 
expression  for  a  point  explosive  source  derived  by  Harkrider  (1964).  It 
can  be  written  as: 


EXPr(w) 


i 

AR(o>,h) 


where: 


EXP  (to) 
R 

/Mb) 

eo 

V(eo) 


AR(w,h) 

r 


Spectral  Amplitude  of  Rayleigh  Wave 

Rigidity  at  depch 

Ellipticity 

Reduced  displacement  potential  in  the  frequency 

W  5/4 

domain;  i^(oo)  =  ty^oo)  (~^)  and  W  is  yield  in 
kilotons 

Wave  number  of  Rayleigh  wave 
Response  of  layered  half  space 
Distance  (km) 


Listed  in  Table  III—  l  are  three  cases  where  a  double  couple 
source  was  fitted  using  three  different  azimuthal  configurations.  The  first 
two  configurations  are  used  in  Sections  IV  and  VI  in  actual  cases,  while  the 
third  could  represent  a  possible  future  case.  The,  optimal  solution  for  each 
case  is  given  in  Table  III- 1,  with  residual  distributions  of  the  source  para¬ 
meters  (see  Appendix  E  for  explanation)  shown  in  Figures  III- 1,  111. 2.  and 
III- 3  for  explosions  l,  2,  and  3  respectively.  If  no  prior  knowledge  of  this 

event  existed,  we  could  assume  an  earthquake  occurred  at  15  km,  with 
.  o  .  o 

dip  of  60  ,  and  slip  of  -30  ,  except  for  the  following: 

•  Examining  the  residual  distributions,  we  see  for  cases  2 

and  3  that  optimal  solutions  exist  at  almost  any  strike 
direction,  indicating  a  highly  circular  pattern.  This  is  not 
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TABLE  IIl-l 

FITTING  THEORETICAL  EXPLOSION  WITH 
DOUBLE-COUPLED  EARTHQUAKE  MODEL 


Case  I.  D. 

Optimal  Solution 

THEO/EXPLOSION/Ol 

Depth 
h  km 

Dip  Angle 
8° 

Slip  Angle 
'  \° 

Strike 

N0°E 

Moment 
1025  dyne-cm 

15 

60 

0 

160 

0.  22  E00 

THEO/ EXPLOSION/ 02 

15 

60 

-30 

20 

0.21E00 

THEO/EXPLOSION/03 

15 

60 

-30 

100 

0.22  E00 

TH  EO  /EX  PLOSIO  N  /  0 1  : 
THEO/ EXPLOSION/ 02: 

THEO  /  EX  xJLOSION  /  0  3 : 


4- station- fit;  4  stations  of  LX/KIRSI/059 
12 -station- fit;  12  stations  of  Mitchell's  Southeastern 
Missouri  Earthquake 
ll- station- fit,  11  VLPE  stations 
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FIGURE  III-l 

RESIDUAL  DISTRIBUTION  FOR  THEORETICAL  EXPLOSION  l 
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true  with  case  1,  where  only  four  azimuths  are  used.  With 
comprehensive  az imuthal  coverage,  though,  complete  indeter¬ 
minancy  in  the  strike  direction  could  indicate  a  source  mis¬ 
match. 

Only  Rayleigh  waves  have  been  used  in  the  fitting  procedure. 
Using  these  source  parameters,  we  generated  both  Rayleigh 
and  Love  wave  radiation  patterns  for  a  double  couple  source 
(Figures  G-lOa,  b  to  G-14a,  b)  from  10  to  50  seconds.  Exam¬ 
ining  the  Rayleigh  wave  patterns,  we  see  a  progressive  elim¬ 
ination  of  lobal  configuration  from  50  to  10  seconds,  with  the 
latter  pattern  being  almost  circular.  At  the  same  time,  the 
Love  wave  patterns  have  sharply  defined  nodes. 

The  circularity  of  the  Rayleigh  wave  radiation  pattern,  though, 
is  not  enough  to  achieve  a  mis-identification;  a  match  must  also 
occur  between  the  spectral  shapes  of  the  explosion  and  theo¬ 
retical  double  couple.  If  we  compare  the  spectra  for  an  explo¬ 
sion  at  shallow  depths  (Figure  III-4)  with  the  Rayleigh  wave 
spectra  (Figure  G-8a)  which  is  approximately  the  same  as  our 
solution,  we  see  that  the  shape  is  similar.  In  Figure  III-5,  a 
comparison  of  the  exact  spectral  shapes  of  these  two  sources 
is  shown  in  detail.  The  spectra  from  an  explosion  at  or  near 
the  surface  essentially  averages  the  variation  of  the  radiation 
patterns  over  the  periods  of  interest.  For  an  explosion  at  4  km, 
this  spectral  averaging  is  not  as  good,  and  it  becomes  worse  for 
deeper  explosions.  The  uniqueness  of  this  case  becomes  evident 
when  one  compares  the  spectral  shapes  of  other  double  couple 
configurations  where  circular  radiation  patterns  exist  with 
shallow  depth  explosion  spectra;  the  shapes  are  found  to  be 
quite  different. 
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Relative  Amplitude  Spectra 


FIGURE  III -5 

COMPARISON  OF  EXPLOSIVE  SOURCE  SPECTRA  AND  A 
PARTICULAR  DOUBLE  COUPLE  EARTHQUAKE 
SOURCE  SPECTRA  (h=15  km,  6=60°,  A  =  -30°) 


Therefore,  from  the  above  points,  we  see  that  misidentification 
occurs  only  in  a  particular  exceptional  case,  and  with  good  azimuthal  coverage, 
this  situation  can  be  clarified.  In  future  studies,  we  will  examine  the  problem 
of  fitting  an  explosive  source  with  an  associated  seismic  relaxation  using  the 
same  double  couple  source. 


SECTION  IV 


I 

*** 


ANALYSIS  OF  THE  SOUTHEASTERN  MISSOURI 
EARTHQUAKE  OF  OCTOBER  21,  1965 

Because  of  the  availability  of  high  quality  Rayleigh  wave  data 
for  the  Southeastern  Missouri  Earthquake  of  October  21,  1965  (m  =5,  2)  a 
good  opportunity  existed  for  testing  the  methodology  under  highly  controlled 
circumstances.  In  this  way,  we  can  determine  what  is  possible  when  we 
analyze  Eurasian  events  if  we  have  extensive  azimuthal  coverage.  The 
Rayleigh  wave  amplitude  spectra  for  this  event  was  obtained  from  Brian 
Mitchell  (personal  communication  -  1973),  having  been  already  corrected  for 
multipathing,  higher  modes,  geometrical  spreading,  instrument  response, 
and  attenuation.  Mitchell's  attenuation  correction  is  comparable  with 
Tryggvason' s  (1965)  for  periods  between  10  and  60  seconds  (note  discussion 
in  Section  VI).  From  the  data,  the  amplitude  spectra  recorded  at  19  obser¬ 
vation  stations  were  chosen  for  the  source  parameter  study.  The  locations 
of  these  stations  with  respect  to  the  epicenter  are  shown  in  Figure  IV- 1. 

Using  the  modified  Tsai'B  method  (Appendix  E)  for  residual 
distributions  and  the  non-linear  least  squares  method  (Appendix  F)  to  deter¬ 
mine  the  absolute  minimum  residual  solution,  parameter  estimates  were 
obtained  for  several  groupings  of  the  stations.  These  solutions  were  obtained 
using  a  double  couple  point  source  placed  in  a  layered  half  space  which  models 
the  central  United  States  (McEvilly,  1963).  A  velocity-density  diagram  of  this 
structure  is  given  in  Figure  D-3,  with  its  Rayleigh  and  Love  wave  group 
velocities  shown  in  Figures  D-5  and  D-6.  The  minimum  residual  statistics 
for  the  12  stations  east  of  the  event  are  shown  in  Figures  IV- 2,  with  the 
solution  deduced  from  the  statistics  given  in  Table  IV- 1.  The  optimal 
solution  is  (with  Mitchell* c  solution  in  parenthesis): 
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VARIATION  OF  e  WITH  RESPECT  TO  DIP  ANGLE 
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VARIATION  OF  c  WITH  RESPECT  TO  SLIP  ANGLE 
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Source  Parameter 


Probable  Region 


%  Confidence 


Focal  Depth:  h  km 

■ 

N 

85 

Dip  Angle:  5° 

40-50 

59.  5 

.  o 

Slip  Angle:  A 

±90 

58 

Strike:  N^°E 

60°- 90° 

75 

Focal  depth  =  4  km  (4  km) 

Dip  Angle  =  40°  or  50°  (50°) 

Slip  Angle  =  ±90°  (+85°  or  +95°) 

Strike  =  N70°E  (N  70°  E) 

25  25 

Seismic  Moment  =  0. 18  x  10  dyne  •  cm  (0. 17  x  10  dyne  •  cm) 

We  see  th?t  the  solution  obtained  from  modified  Tsai's  method 
and  Mitchell's  solution  are  in  very  close  agreement.  The  distribution  of 
minimum  residuals  gives  us  some  idea  of  confidence  regions  for  our  parameter 
choices  (Figures  IV-2).  For  example,  from  Figure  IV-2a,  examination  of  the  vari¬ 
ation  of  the  minimum  residual  with  respect  to  depth  yields  the  following  results: 

•  Fixing  the  dip,  slip,  and  strike  at  the  optimal  solution, 
variation  over  depth  yields  a  minimum  residual  at  4  km. 

•  Varying  all  of  the  parameters,  the  minimum  of  the  minimum 
residuals  at  each  depth  occurs  at  4  km.  A  standard  minimum 
residual  is  then  chosen  to  determine  the  frequency  solutions  that 
occur  at  certain  depths  with  minimum  residuals  below  the 
standard  value. 

•  As  displayed  in  Figure  IV-2a  and  listed  in  Table  IV-l,  85 
per  cent  of  the  models  with  minimum  residuals  less  than  the 
standard  have  depths  between  two  and  four  kilometers. 

As  mentioned  earlier,  solutions  were  obtained  for  several 
combinations  of  stations.  Figures  IV- 3  display  spectra  fit  at  six  stations 
chosen  at  random.  We  see  that  the  19-station  fit  is  almost  identical  to  the 
12 -station  fit.  Because  the  amplitude  of  the  theoretical  double  couple  is  sym¬ 
metrical  about  180°,  the  additional  seven  stations  are  not  needed.  Good 
coverage,  such  as  this,  in  two  successive  quadrants,  produce  excellent 
results.  Successive  solutions  were  obtained  by  reducing  the  number  of 
stations.  The  solution  which  used  the  fewest  number  of  stations,  and  still 
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retained  the  same  depth  estimate,  contained  only  one  station  in  each  of  the 

eastern  quadrants.  These  results  indicate  what  is  possible  with  high  quality 
data. 

Application  of  the  non-linear  regression  analysis  yielded  the 
following  solution; 


Focal  depth  «4,0  km 
F»ip  Angle  =  53.  1° 

Slip  Angle  -  -79.  7° 

Strike  -  N80.  4°E 

O  C 

Seismic  Moment  =  0.  167  x  10  dyne  •  cm 

(mb  =  5.4  using  w2-model  [Tsai,  19721). 

This  solution  represents  the  absolute  minimum  residual  fit. 

A  radiation  pattern  plot  for  periods  10-50  seconds  (AT  =  5  seconds)  is  shown 
in  Figures  IV-4.  The  fits  at  T=25,  30,  and  35  seconds  are  especially  close. 

The  consistency  of  the  depth  estimate  using  the  previous 
techniques  is  quite  encouraging.  Mitchell  (1973),  using  the  same  data  in 
conjunction  with  a  trial-and-error  technique  guided  by  partial  derivatives  of 
the  source  parameters,  obtained  very  much  the  same  result.  Using  body 
wave  data,  Kisslinger  and  Nultti  (1965)  obtained  depths  of  <6  km. 


Before  we  proceed  to  data  with  less  control,  we  should  examine 
the  effect  that  the  layered-half  space  velocity  structure  has  on  the  parameter 
estimates.  As  stated  earlier,  the  solutions  for  this  event  have  been  obtained 
using  the  Central  United  States  model.  Three  other  plane  layered  models, 
namely  Gutenberg  .  Bullen,  Tibet  Plateau,  and  Hamilton-Healy,  were  used  in 
processing  the  same  data.  The  velocity  structure  and  group  velocity  curves 
for  these  models  is  given  in  Appendix  D.  All  four  earth  models  are  for 
continental  areas,  with  the  main  differences  in  the  upper  crust.  Source 
parameter  estimates  using  these  four  models  are  listed  in  Table  IV-2. 
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FIGURE  IV-4a 
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FIGURE  IV-4b 

RADIATION  PATTERNS  OBTAINED  USING  NON-LINEAR  REGRESSION-- 
THE  SOUTHEASTERN  MISSOURI  EARTHQUAKE  OF  OCTOBER  21,  1965 
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THE  SOTTTHEASTERN  MISSOURI  EARTHQUAKE  OF  OCTOBER  21,  1965 
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RADIATION  PATTERNS  OBTAINED  USING  NON-LINEAR  REGRESSION-- 
THE  SOUTHEASTERN  MISSOURI  EARTHQUAKE  OF  OCTOBER  21,  1965 
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RADIATION  PATTERNS  OBTAINED  USING  NON-LINEAR  REGRESSION- - 
THE  SOUTHEASTERN  MISSOURI  EARTHQUAKE  OF  OCTOBER  21,  1965 
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TABLE  IV- 2 

SOUTHEASTERN  MISSOURI  EARTHQUAKE 
VARIATION  OF  SEISMIC  SOURCE  PARAMETER 
ESTIMATION  WITH  RESPECT  TO  DIFFERENT  EARTH  MODELS 


Earth  Model 

Source  Parameter  Estimation 

Focal  Depth 
h  km 

Dip  Angle 
8° 

Slip  Angle 

Strike 

N<f>°E 

Moment 

M  dy  cm 

Central  U.  S. 

4 

40 

50 

±90° 

70 

G.  18 

56km-Tibet  Plateau 

3 

40 

-60° 

55 

0. 12 

Gutenberg 

4 

40 

50 

*90° 

no 

0.18 

Hamilton- Healy 

5 

40 

50 

±90° 

70 

0.13 
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From  these  results,  one  comes  to  the  conclusion  that  structure 
has  little  effect  on  the  depth  estimate.  This  result  is  quite  puzzling.  The 
large  differences  between  the  Tibet  Plateau  high  velocity  crust  and  the 
Hamilton- Healy  sedimentary  crust  yield  only  a  2  km  change  in  the  depth 
estimate.  Future  studies  will  examine  the  effect  layer  parameters  have  on 

depth  estimation. 


IV- 2  3 


SECTION  V 

ANALYSIS  OF  TWO,  LOW  MAGNITUDE  ITALIAN  EVENTS 


In  the  analysis  of  the  Southeastern  Missouri  Earthquake,  we 
were  able  to  obtain  a  depth  estimate,  using  only  two  stations  in  consecutive 
quadrants,  which  was  consistent  with  the  depth  estimate  obtained  using  12 
stations.  Since  this  situation  (a  paucity  of  stations)  usually  is  present  in  the 
analysis  of  Eurasian  data,  examination  of  just  such  a  problem  was  under¬ 
taken  where  some  control  of  the  source  parameters  was  available. 

Two  low'  magnitude  (m^=4.1)  events  occurred  off  the  east 
coast  of  Italy  (Figure  V-l),  approximately  one  day  apart  (Table  V-2),  and 
were  recorded  by  three  VLPE  stations  (Table  V-3).  Our  control  of  the  source 
parameters  lie  in  the  body  wave  solution  of  an  event  in  the  same  tectonic 
region  (McKenzie,  1972)  along  the  Italian  east  coast  (Figure  V-2).  The  crust 
structure  for  the  source  region,  obtained  from  Closs  (1969),  is  shown  in 
Table  V-l. 

Corrections  from  the  source  to  the  three  common  VLPE  sites 
were  applied  using  Tsai's  attenuation  correction  (1972)  for  the  all  paths;  to 
TLO,  OGD,  and  KON.  A  series  of  narrow  band  filters  were  applied  to  the 
waveforms  to  determine  if  multipathing  was  present  (see  Appendix  B  for 
method),  and  very  little  was  found. 

Because  of  the  uneven  quality  of  the  data,  we  only  applied  the 
modified  Tsai's  method,  with  the  accuracy  furnished  by  the  residual  distri¬ 
butions  thought  to  be  sufficient.  The  minimum  residual  statistics  are  shown 
in  Figures  V-*  for  CITLY/134,  and  in  Figures  V-4  for  CITLY/143  using  all 
three  VLPE  stations.  A  summary  of  these  statistics  is  given  in  Table  V-4. 


LOCATIONS  OF  TWO,  LOW  MAGNITUDE  ITALIAN  EVENTS 
RECORDED  AT  THREE  VLPE  STATIONS 
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TABLE  V- 3 


OBSERVATION  STATIONS  FOR  LX/CITLY/ 


1.  Event:  LX/CITLY/143 


Site  I.  D. 

Azimuth  from  Source 
(+  for  N-E) 

A  in  Degree 

A  in  Km 

4  TLO 

-100. 3 

13.  5 

1500.  3 

6  KON 

-7.  0 

16.  1 

1788. 7 

7  OGD 

-58.  7 

61.  8 

6858. 5 

2.  Event:  LX/CITLY/134 


Site  I.  D. 

Azimuth  from  Source 
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A  in  Km 
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VARIATION  OF  €  WITH  RESPECT  TO  DEPTH 


Dip  Angle 


•It 


%  Freq 


LX/CITLY/143 


20.0 


>. 

o 

c 

4) 
0 
c r 

4) 


10.  0 


FIGURE  V-4b 

VARIATION  OF  6  WITH  RESPECT  TO  DIP 


V-ll 


I 


TABLE  V-4 

SOLUTIONS  FOR  ITALIAN  EVENTS 


Analyzed  Individually 

_LX/CITLY  /  1  34  (2/4/72)1  .LX/CITLY/143  (2/5/72) 


Parameters 
h  (km) 

8 

A 

<t> 


Value 

Statistics 

Value 

Statistics 

22 

20-26  ;  88% 

20 

18-24  ;  81% 

50° 

40-60°;  68% 

50° 

50°-  70°;  65% 

-30° 

— 6(^°  30°  89% 

0° 

- 30 °- 30°;  99% 

N110°E 

21 

1.2x10 

80°- 110°;  34% 

N105°E 

1.4xl021 

9 5°- 115°;  3  7% 

I 


Analyzed  Using  Spectral  Ratio's  of  Events 


T 

S  «j) 


Trial  Event 


T 

S  (U) 


Ref.  Event 


ROBM, 

R°B((j), 


mTE  Cf  :TE  STK(I 

!Mre  OM  resCB(c 

RE 


LX/CITLY/134  With  LX/CITLY/143  As  Referenc< 


Parameters 

h  (km) 

5 

A 

* 


Value 

18 


N90°E 


Statistics 

16-20  ;  87% 
40  -70  ;  74% 
±30°;  61% 
80°-ll0°;  34% 


M  (dyne-cm) 


1.  3x10 


The  source  parameters  for  each  event  were  similar,  which 
should  be  expected  because  of  their  proximity  (<l°).  The  residual  distri¬ 
bution  for  depth  was  especially  encouraging,  with  80  per  cent  of  the  solutions 
lying  within  a  6  km  range  for  both  events.  The  dip  and  strike  for  both  events 
are  in  close  agreement  with  McKenzie's  (1972)  solution  for  the  event  in  the 
same  tectonic  region.  The  seismic  moments  also  correlate  with  the  body 

wave  magnitudes  for  each  event,  yielding  m  =  3.8  to  4.  0  using  Tsai's 

2  b 

co  -model  (1972). 

The  almost  identical  location  of  these  events  offered  a  good 
opportunity  to  determine  eource  parameters  using  spectral  ratios.  It  is 
hoped  that  this  method  will  enable  us  to  examine  low  magnitude  events  using 
large  magnitude  events  from  the  same  region  as  a  reference  source.  The 
spectral  ratio  is  based  on  the  assumption  that  one  or  more  events  from  a 
region  of  interest  have  been  rigorously  analyzed,  with  the  source  parameter 
estimates  obtained  and  corroborated  from  independent  sources.  Then  small 
events  in  the  same  region,  recorded  at  the  same  stations  so  source-station 
travel  paths  are  similar,  can  be  determined  by  essentially  comparing  source 
excitation  functions.  Examining  the  equation  in  Table  V-4,  we  see  that  the 
path  transfer  functions  (G(w))  and  instrument  response  (I (oj) )  will  cancel  if 
we  have  the  same  source- station  paths.  When  CITLY/134  was  analyzed  using 
CITLY/143  as  a  reference,  we  obtained  parameter  estimates  reasonably 
close  to  those  determined  using  only  the  event  itself  (Table  V-4).  The  depth 
estimate  changed  4  km. 

The  spectral  ratio  method  can  also  give  us  some  indication  of 
the  effectiveness  of  our  path  correction  if,  as  in  this  case,  individual  esti¬ 
mates  of  the  source  parameters  are  known.  Since  multipathing  was  not 
observed,  our  correction  for  effective  Q  must  be  inaccurate  for  one  or  more 
paths.  This  is  indicated  if  we  examine  the  spectral  fits  obtained  using  the 
individual  event  parameter  estimates.  Figures  V-5  show  fits  for  CITLY/134, 


Amplitude  (  Cm-Sec  )  Amplitude  (  Cm-Sec 


and  Figures  V-6  show  the  same  for  CITLY/143.  For  both  events,  we  see 
that  the  spectral  match  at  TLO  breaks  down  completely  at  longer  periods. 
Evidently,  our  effective  Q  is  not  modeling  the  actual  attenuation.  Future 
studies  will  be  directed  at  obtaining  an  accurate  estimate. 

Finally,  two  other  cases  of  interest  were  analyzed,  and  the 
results  are  summarized  below: 

•  Using  Event  CITLY/134  as  a  reference,  CITLY/143  was 
analyzed  using  spectral  ratios.  The  depth  estimate  was  22  km, 
with  the  possible  range  22  to  26  km  containing  80  per  cent 

of  the  models. 

•  Both  events  were  analyzed  individually  using  only  two  stations, 
KON  and  OGD.  The  depth  estimate  for  both  events  was  24  km, 
but  the  other  parameters  varied  a  great  deal  from  the  three 
station  fit,  with  the  seismic  moments  for  both  events  about 

a  third  larger. 
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FIGURE  V-6a 

SPECTRAL  AMPLITUDE:  LX/CITLY/ 143:  KON 
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FIGURE  V-6b 

SPECTRAL  AMPLITUDE:  LX/CITLY/ 143:  OGD 
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FIGURE  V-6c 

SPECTRAL  AMPLITUDE:  LX/CITLY/143:  TLO 


SECTION  VI 


EARTHQUAKES  IN  SINKIANC 


Having  analyzed  events  with  varying  degrees  of  control  on 
the  source  parameters,  we  will  now  analyze  six  earthquakes  in  Sinkiang, 
China,  for  which  we  have  no  reliable  independent  estimates  of  the  source 
parameters.  These  events  have  been  chosen  from  about  thirty  events  which 
originated  in  that  region  during  1972,  the  sole  criteria  being  that  they  were 
observed  at  two  or  more  VLPE  stations.  The  locations  of  these  events  are 
displayed  in  Figure  VI-1,  with  their  locations,  date  and  origin  time,  and 
body  wave  magnitude  given  in  Table  VI-1.  Information  concerning  the  obser¬ 
vation  stations  is  given  in  Table  VI-2  for  each  event. 

This  section  is  divided  into  three  parts.  In  the  first  we  will 
discuss  attenuation  and  multipathing  corrections.  The  second  part  contains 
solutions  and  residual  distributions  for  each  event,  and  the  third  part  presents 
solutions  obtained  through  spectral  ratios. 

A.  PATH  CORRECTIONS 

Narrow  band  filter  analysis,  as  discussed  in  Appendix  B,  was 
applied  to  every  Rayleigh  and  Love  wave  record.  Examples  of  the  narrow 
band  output  are  given  in  Appendix  H,  with  the  associated  group  velocity 
curves.  Observed  Rayleigh  wave  amplitude  spectra  was  available  for  all 
eight  sampled  periods  (T  =  15  to  50  seconds,  4T  =  5  sec)  of  the  six  events 
at  most  of  the  observation  stations;  while  for  the  Love  waves,  the  observed 
amplitude  spectra  was  incomplete  at  most  stations.  Except  for  event 
LX/KIRSI./059,  the  spectral  amplitudes  are  relatively  small.  For  station 
CHG,  though,  the  spectral  amplitudes  are  unreasonably  small,  indicating 
unusual  attenuation.  (The  source- station  path  of  Sinkiang  to  CHG  is  shown 
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TABLE  VI-l 

SINKIANG  EARTHQUAKES 


T 

i 


Fvent  I.  D. 

Location 

Date  and  Time 

Latitude 

Longitude 

LX/KIRSI/05  9 

40.  3N 

79.  0E 

01/15/72 

20.  21.  50 

LX/SSINK/  990 

39.  4N 

75.  3E 

12/03/72 

08.  54.47 

LX/NSINK/ 008 

41.  8N 

84.  5E 

01/02/72 

10.  27.  35 

LX/SINKI/178 

41.  7N 

80.  7E 

02/16/72 

23.  19.20 

LX/NSINK/566 

44.  6N 

81.  IE 

07/05/72 

01.  09.  53 

LX/SINKI/164 


3  9.  9N 


77.  4E 


02/11/72 
05.  55.46 


TABLE  VI-2 

OBSERVATION  STATIONS  FOR  SIUKIANG  EARTHQUAKES 


Fvent  I.  D. 

Site 

I.  D. 

Azimuth  To  Source 

A  In  Degree 

A  In  km 

(+  for  N-E) 

LX/KIRSl/059 

2 

CHG 

135.4 

27.  5 

3056.  3 

3 

FBK 

19.  5 

68.  7 

7648. 4 

4 

TLO 

-60. -6 

61.0 

6790. 0 

7 

OGD 

-19.  7 

95.  3 

10604.  7 

LX/SSINK/990 

2 

CHG 

128.2 

29.0 

3226.  0 

6 

KON 

-40.  7 

45.  0 

5004.  0 

9 

ALQ 

1.5 

105.  9 

11775.  9 

LX/NSINK/  008 

2 

CHG 

147.  3 

26.  1 

2896. 9 

4 

TLO 

-58.  8 

63.  8 

7105.  3 

LX/SINKI/178 

3 

FBK 

20.  3 

67.  0 

7453. 5 

6 

KON 

-41.  7 

45.  9 

5113.  6 

LX/NSINK/566 

4 

TLO 

-61.  7 

mm 

67U.  0 

8 

KIP 

53.  2 

MSM 

10563.  4 

LX/SINKI/164 


3  FBK 
6  KON 


18.  9 
-40.  9 


69.  5 
45.  6 


7735.  I 
5079.  8 


in  Figure  VI-2,  as  well  as  paths  for  all  the  VLPE  stations.)  This  high 
attenuation  also  shows  up  in  the  amplitude  spectra  for  explosions  in  Eastern 
Kazakh.  Listed  in  Table  VI- 3  are  event  information  for  two  of  these. 

The  explosions  were  used  in  an  attempt  to  obtain  an  effective 
Q  for  the  Sinkiang  -  CHG  travel  path  from  Tryggvason's  method  (Appendix  A). 
His  method  assumes  that  equal  energy  is  transmitted  in  all  directions  from 
the  source,  which  is  ideally  suited  for  explosion  data.  In  Figure  VI-3  are 
shown  energy  attenuation  coefficients  kg  for  the  two  explosions  determined  by 
Tryggvason's  method.  The  results  are  hard  to  interpret,  mainly  because 
his  method  needs  good  azimuthal  coverage  and  purely  continental  paths.  We 
tried  to  adjust  the  amplitudes  of  EKAZH/013  using  Von  Seggern's  (1971) 
amplitude  equalization  technique  for  estimating  differences  between  contin¬ 
ental  and  oceanic  paths,  and  then  apply  Tryggvason's  technique.  Again,  the 
results  are  difficult  to  interpret  for  the  reasons  given  above.  We  finally 
decided  to  implement  a  two- station  correction  (Appendix  A)  between  CHG 
and  ALQ  using  Tryggvason's  energy  attenuation  coefficient.  Although  CHG 
and  ALQ  are  not  on  the  same  great  circle  path  through  the  source  area,  we 
are  using  explosion  data,  and  the  ALQ  travel  path  is  mostly  continental  (as 
shown  consistently  by  group  velocity  curves  from  the  Sinkiang  region).  The 
results  of  this  technique  are  shown  in  Figure  VI-4.  The  energy  attenuation 
coefficient  for  CHG  is  quite  high,  indicating  high  absorption  over  the 
frequencies  of  interest. 

B.  AMPLITHDE  SPECTRAL  FITTING  FOR  INDIVIDUAL  EVENTS 

As  with  the  Italian  data,  the  quality  is  highly  variable,  and 
we  only  applied  the  modified  Tsai's  method,  with  the  accuracy  furnished  by 
the  residual  distributions  thought  to  be  sufficient.  The  optimal  solutions  of 
source  parameters  for  these  six  earthquakes  are  listed  in  Table  Vl-4,  and 
the  probable  ranges  with  percentage  confidences  of  the  source  parameters 
are  listed  in  Table  VI-5.  The  distributions  of  minimum  residuals  oi  each 
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SOURCE-STATION  PATHS  FOR  THE  VLPE  STATIONS  FROM  SINKIANG  EVENTS 


TABLE  VI- 3 

TWO  EXPLOSIONS  FROM  EASTERN  KAZAKH 


LX/EKAZH/013:  Location:  49.  8  N  Lat.  ,  78. 1  E  Long. 

Magnitude:  =  5.7 

Date:  12/10/72 

Origin  Time:  04:26:58 

Recording 

Station 

Location 

Azimuth 
From  Source 

A  (km) 

Lat. 

Long. 

CHG 

18.  8  N 

99.  0  E 

144.  0° 

3911.  2 

KIP 

21.  4  N 

15b.  0  W 

50.  6° 

10385.  8 

ALQ 

34.  9  N 

106.  5  W 

3.  8° 

10611.  3 

ZLP 

16.  5  S 

68. 1W 

-51.1° 

15233.3 

LX/EKAZH/876:  Location:  49.  9N  Lat.  ,  78.  8  E  Long. 


Magnitude:  =  6.  2 

Date:  11/2/72 
Origin  Time:  01:26:58 


Recording 

Station 

Location 

Azimuth 
From  Source 

A  (km) 

Lat. 

Long. 

KON 

59.  7  N 

9.  6  E 

-'t 

oo 

-'t 

i 

4364.  0 

KIP 

21.  4N 

158.  0  W 

51.2 

10339.  7 

ALQ 

34.  9  N 

106.  5  W 

4.4 

10596.  7 

Tryggvaaon's  Fnergy  Attenuation  Coefficient  k 


ENERGY  ATTENUATION  COEFFICIENT  FOR  CHG 


TABLE  VI-4 

ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  AMPLITUDE  SPECTRAL  FITTING 
BASED  ON  MINIMUM-RESIDUAL  CRITERION 


Event  I.  D. 

Optimal  Solution 

Depth 

h  km 

Dip  Angle 

5  Degree 

Slip  Angle 

A  Degree 

Strike 

N0°E 

Moment 

25  . 

10  dyne -cm 

LX/KIRSI/059 

42.  5 

40.  0 

0.  0 

170.  0 

0.  768E-01 

LX/SSINK/990 

47.  5 

50.0 

-60.  0 

120.  0 

0. 103E-02 

LX/NSINK/008 

4.0 

30.  0 

30.  0 

165.  0 

0.  669E-03 

LX/SINKI/178 

10.  0 

40.  0 

60.  0 

150.  0 

0.  452E-03 

LX/NSINK/566 

14.  0 

60.  0 

30.  0 

20.  0 

0.  338E-03 

LX/SINKI/164 

6.  0 

40.  0 

30.  0 

150.  0 

0.  512E-03 
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TABLE  VI- 5 

ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  AMPLITUDE  SPECTRAL  FITTING 
BASED  ON  DISTRIBUTION  OF  MINIMUM  RESIDUAL  CRITERION 


Event  I,  D. 

Source  Parameters  ; 

Depth  h  km 

Dip  Angle  8° 

Slip  Angle  A° 

Strike  N$°E 

Probable 

_  1 
Range 

% 

Confi¬ 

dence 

Probable 

_  2 
Range 

% 

Confi¬ 

dence 

Probable 

„  3 

Range 

% 

Confi¬ 

dence 

Probable 

4 

Range 

% 

Confi¬ 

dence 

30-40 

46 

-60 

25 

75-100 

37 

LX/KIRSI/059 

40-52.5 

83 

70-90 

37 

0-60 

55 

160-180 

25 

115-125 

23 

LX/SSINK/990 

42. 5-55 

78 

40-60 

90 

30,  60 

71 

155-180 

41 

30 

30 

ESI 

42 

LX  /N  SINK/  008 

4-8 

85 

60-80 

55 

fiSS 

48 

100-135 

48 

LX/SINKI/178 

6-12 

80 

30-60 

80 

30,  60 

83 

150-170 

20 

15-35 

mm 

LX/NSINK/566 

8-18 

75 

60-80 

75 

30 

78 

105-125 

S9 

30-40 

54 

LX/SINKI/164 

2-8 

89 

70-80 

_ 

25 

30,  60 

75 

150-165 

39 

1.  Ah  =  2  km,  0<h<20  km,  Ah  s  2.  5  km,  h>20  km 

2.  AS  =  10° 

3.  A  K  =  30° 

4.  A  <f>  =  5° 
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parameter  of  each  event  is  shown  in  Figures  VI-5  to  VI-10  with  their 
spectral  fits. 

Referring  to  Table  VI-4,  these  six  events  may  be  divided 

into  two  groups  according  to  their  focal  depth;  (l)  LX/KIRSI/05  9  and 

LX/SSINK/990  have  the  deeper  focal  depth  at  about  45  km;  (2)  the  other  four 

have  shallow  focal  depths  around  5  to  15  km.  From  the  dip  angle  estimation, 

these  six  earthquakes  have  fault  planes  dipping  from  30°  to  60°.  The  strikes 

of  all  but  one  event,  LX/NSINK/566,  are  oriented  approximately  N  150°  E. 

The  magnitudes  of  the  seismic  moments  estimated  for  these  events  agree 

with  their  m,  values  from  the  PDE  list  in  the  sense  that  the  larger  the  m 
b  b 

values  are,  the  larger  the  seismic  moments  will  be  for  the  same  focal  depth. 

In  general,  the  percentage  confidence  of  the  focal  depth 
estimation  is  about  80%  for  about  6  km  uncertainty  in  the  shallow  focal 
depths  and  about  10  km  uncertainty  in  the  deeper  focal  depths.  The  resolution 
of  the  rest  of  the  source  parameters  is  mixed,  with  the  strike  distributions 
generally  poor.  This  might  be  due  to  the  insufficient  azimuthal  coverage  and 
incomplete  Love  wave  spectra. 

Visually  judging  the  spectral  fits,  those  for  the  four  shallow 
earthquakes  are  reasonably  good  at  both  stations,  while  those  for  the  two 
deeper  events  are  mixed  and,  in  general,  are  just  fair. 

C.  AMPLITUDE  SPECTRAL  FITTING  USING  SPECTRAL  RATIOS 

Among  the  six  Sinkiang  earthquakes,  there  are  two  pairs  of 
events  which  have  common  observation  stations.  Referring  to  Table  VI-2, 
they  are:  (l)  LX/KIRSI/059  and  LX/NSINK/008;  (2)  LX/SINKI/178  and 
LX/SINKI/164,  Common  observation  stations  for  the  first  pair  are  CHG 
and  TLO;  and  those  for  the  second  pair  are  FBK  and  KON.  Both  pairs  of 
events  occur  within  5°. 
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FIGURE  VI- 5a 

RESIDUAL  DISTRIBUTIONS  FOR  LX/KIRSI/059 
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RAYLEIGH  WAVE  SPECTRAL  FITS  FOR  LX/SSINK/990 
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FIGURE  VI-6c 

LOVE  WAVE  SPECTRAL  FITS  FOR  LX/SSINK/990 


VI-18 


FIGURE  VI- 7a 

RESIDUAL  DISTRIBUTIONS  FOR  LX/NSINK/008 
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FIGURE  VI- 7b 
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RAYLEIGH  WAVE  SPECTRAL  FITS  FOR  LX/NSINK/008 
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FIGURE  VI- 7c 

LOVE  WAVE  SPECTRAL  FITS  FOR  LX/NSTNK/008 
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RESIDUAL  DISTRIBUTIONS  FOR  LX/SINKI/178 
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FIGURE  VI-8c 

LOVE  WAVE  SPECTRAL  FITS  FOR  LX/SINKI/178 
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FIGURE  VI- 9a 

RESIDUAL  DISTRIBUTIONS  FOR  LX/NSINK/566 
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FIGURE  VI- 9b 

RAYLEIGH  WAVE  SPECTRAL  FITS  FOR  LX/NSINK/566 
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FIGURE  VI- 9c 

LOVE  WAVE  SPECTRAL  FITS  FOR  LX/NSINK/566 
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FIGURE  VI-10a 

RESIDUAL  DISTRIBUTIONS  FOR  LX/SINKI/164 
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FIGURE  VI-lOc 

LOVE  WAVE  SPECTRAL  FITS  FOR  LX/SINKI/164 


The  estimations  of  source  parameters  for  events  LX/NSINK/008, 
LX/SINKl/178 ,  and  LX/SINKI/164  using  ratio-of-events  fitting  with  minimum 
residual  and  distribution-of-minimum-residual  criteria  are  given  in  Tables 
VI-6,  VI-7,  and  VI-8  respectively.  The  residual  distributions  are  also 
displayed  in  Figures  VI-11,  VI-12,  and  VI-13.  From  these  results,  we  can 
reach  the  following  conclusions: 

•  For  event  LX/NSINK/008,  the  optimal  solution  of  source 

parameters  obtained  from  the  ratio-of-events  fitting  is  some¬ 
what  different  from  that  obtained  with  the  individual  spectral 
fit,  except  for  the  focal  depth  estimation.  This  discrepancy 
may  be  due  to  two  reasons:  (l)  the  locations  of  the  two  events 
are  not  close  enough,  or  (2)  the  original  parameter  estimation 
of  LX/KIRSI/059,  the  reference  event,  may  be  not  good 
enough. 


•  For  events  LX/SINKI/178  and  LX/SINKI/164,  the  source 

parameter  estimations  from  spectral  ratios  are  quite  compar¬ 
able  with  those  from  amplitude  spectral  fitting.  This  seems 
to  firm  the  good  spectral  fits  obtained  for  the  shallow  events. 


TABLE  VI- 6 


LX/NSINK/008 

ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  RATIO- OF- EVENTS  FITTING 

REFERENCE  EVENT:  LX/KIRSI/059 

SOLUTIONS  FOR  EVENT:  LX/NSINK/008 

A.  Solution  by  Minimun-Residual  Criterion 


Optimal  Solution 

Event  I.  D. 

Depth 

Dip  Angle 

Slip  Angle 

Strike 

Moment 

h  km 

5  Degree 

X  Degree 

N0°E 

,„25  J 

10  dyne- cm 

LX/NSINK/008 

90.0 

±90.0 

0. 126E-02 

B.  Solution  by  Distribution-of-Minimum- Residual  Criterion 


Event  I.  D. 

LX/NSINK/008 

Source  Parameter 

Probable  Range 

%  Confidence 

h 

6-10 

86 

5 

70-90 

83 

X 

±90 

74 

0 

125-130 

47 
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FIGURE  VI- II 

RESIDUAL  DISTRIBUTIONS  FOR  LX/NSINK/008 
DETERMINED  BY  RATIO  OF  EVENTS 
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TABLE  VI-7 
LX/SINKI/178 

ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  RATIO- OF- EVENTS  FITTING 


REFERENCE  EVENT:  LX/SINKI/164 


SOLUTIONS  FOR  EVENT:  LX/SINKI/178 

A.  Solution  by  Minimun- Residual  Criterion 


Event  I.  D. 

Optimal  Solution 

Depth 

h  km 

Dip  Angle 

8  Degree 

Slip  Angle 

X  Degree 

Strike 

N  </>°E 

Moment 
lrt25  - 

10  dyne- cm 

LX/SINKI/178 

12.  0 

30,  60 

±90.0 

125.0 

0.  377E-03 

B.  Solution  by  Distribution-of-Minimum- Residual  Criterion 


Event  I.  D. 

LX/SINKI/r8 

Source  Parameter 

Probable  Range 

%  Confidence 

h 

IC-14 

85 

8 
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TABLE  VI-8 


LX/SINKI/164 


~~ i  1MATIONS  OF  SOURCE  PARAMETERS 
TAINED  BY  RATIO-OF-EVENTS  FITTING 


REFERENCE  EVENT:  LX/SINKI/178 

SOLUTIONS  FOR  EVENT:  LX/SINKI/164 

A.  Solution  by  Minimun-Residual  Criterion 


Event  I.  D. 

r - -  - Optimal  Solution 

Depth 

b  km 

Dip  Angle 

5  Degree 

Slip  Angle 

X  Degree 

Strike 

N0°E 

LX/SINKI/164 

— - — — — - _i 

WiM 

60.  0 

30.  0 

B.  Solution  by  Distribution-of-Minimum-Residual  Criterion 


Event  I.  D. 

LX/SINKI/164 

Source  Parameter 

Probable  Range 

"  H 

%  Confidence 

h 

4-8 

78 

6 

50-80 

83 

X 

.130 

70 

<P 

155-160 

57 

LX/SINKI/lfet 


FIGURE  VI-13 

RESIDUAL  DISTRIBUTIONS  FOR  LX/SINKI/164 
DETERMINED  BY  RATIO  OF  EVENTS 
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SECTION  VII 


CONCLUSIONS  AND  DISCUSSIONS 

The  following  remarks  summarize  the  results  in  the  pre- 
ceeding  sections  concerning  source  parameter  estimates  obtained  from 
surface  wave  data  using  spectral  fitting  procedures,  both  exhaustive  and 
iterative: 

•  In  an  effort  to  examine  the  uniqueness  of  our  parameter 
estimates  under  the  minimum  residual  criteria,  initial  studies 
yielded  the  following  results:  (l)  spectral  shape  is  most 
sensitive  to  variations  in  depth;  (2)  No  readily  observable  general 
spectral  property  characterizes  a  source  parameter's  value. 
Holes  in  the  Rayleigh  wave  spectrum  exist  only  for  almost 
vertical  strike  slip  faults;  (3)  A  double  couple  source  applied  to 
theoretical  Rayleigh  wave  data  yields  a  non-shallow  depth  esti¬ 
mate,  out  from  the  distribution  of  possible  strike  angles,  radia¬ 
tion  pattern  shape,  and  spectral  shape,  we  see  that  this  solution 
ir  an  exception.  Future  studies  will  examine  the  source  des¬ 
cription  and  layered  half  space  physical  properties  to  see  what 
parameter  affects  spectial  shape  the  most. 

•  Analysis  of  the  high  quality  Rayleigh  wave  data  from  the  South¬ 
eastern  Missouri  Earthquake  of  October  21,  1965,  yielded 
good  parameter  corroboration  (4  km  depth  vs.  <  6  km  body 
wave  estimate)  and  excellent  spectral  fits.  Although  other 
parameter  estimates  varied,  the  depth  estimate  did  not  change 
for  two  stations  in  successive  quadrants.  Parameter  estimates 
were  obtained  using  four  different  layered  earth  models,  and 


yielded  the  puzzling  result  of  very  little  effect  on  the  depth 
estimate  (Ah  =  2  km),  although  the  other  parameters,  including 
seismic  rrorr  ant,  varied  a  great  deal.  Future  studies  will 
examine  this  problem  in  detail. 

The  analysis  of  two,  low  magnitude  4. 1)  Italian  events 

yielded  good  corroboration  with  body  wave  data.  These 
results  were  obtained  using  three  VLPE  stations.  Each  surface 
wave  record  was  analyzed  for  multipathing  u.iing  narrow  band 
filters,  and  very  little  of  it  was  found.  This  method  has 
proved  to  be  successful,  and  the  results  are  easy  to  interpret. 
These  events  were  also  analyzed  using  spectral  ratios.  This 
method  assumes  that  master  events,  with  accurate  source 
parameter  estimates,  have  been  obtained  for  the  region  of 
interest.  Therefore,  by  taking  spectral  ratios,  travel  path 
effects  have  been  eliminated.  Comparison  of  this  technique 
with  the  individual  fits  were  favorable.  One  problem,  though, 
arose  in  the  spectral  fit  at  TLO.  An  attenuation  study  of  this 
travel  path  will  be  undertaken. 

Six  events  from  Sinkiang,  recorded  at  VLPE  stations  during 
1972,  were  analyzed.  Because  of  the  mixed  quality  of  the  data, 
the  minimum  residual  distribution  statistics  were  particularly 
helpful  in  analyzing  the  optimal  solution.  Severe  multipathing 
rendered  some  Love  wave  data  useless.  The  effective  Q 
corrections  are  thought  to  be  inadequate,  especially  for  the 
CHG  travel  path.  The  source  parameter  analysis  yielded 
mixed  results,  with  the  four  shallow  events  thought  to  have  the 
better  solutions.  Greater  azimuthal  coverage  should  reduce 
much  of  the  uncertainty.  Future  studies  will  include  the 
analysis  of  seismic  regions  within  Eurasia  using  VLPE  data. 
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APPENDIX  A 

TWO  METHODS  FOR  DETERMINING  AN  EFFECTIVE  Q 


All  methods  for  correcting  for  effective  attenuation  involves  a 
correction  for  geometrical  spreading,  which  can  be  written: 


A  =  A 

/Re  sin  Ap  \ 

crQ  op 

\  To  1 

where:  A 

cro 

-  Amplitude  corrected  to  field  point  r 

o 

A 

op 

-  Observed  amplitude  at  field  point  p 

ap 

-  Distance  between  source  and  field  point  p 

R 

e 

-  Radius  of  Earth 

The  geometrical  spreading  correction  is  included  in  the  equations  of  both  of 
the  following  corrections. 

TRYGGVASON'S  APPROACH 

The  basic  assumption  in  Tryggvason's  (1965)  approach  is  that 
equal  energy  is  transmitted  in  all  directions  from  the  source.  Therefore,  this 
method  is  particularly  useful  when  applied  to  explosion  data.  It  is  also 
assumed  that  wave  energy  is  not  transferred  from  one  frequency  band  to 
another.  In  other  words,  spectral  modulations  resulting  from  multipathing 
have  been  eliminated.  The  total  energy  in  a  narrow  frequency  band  at  distance 
D  from  the  source  is: 


E 

o 


-nk  z  .  . 

e  =  C. A  sin  A 
>  l 


(A-Z) 


A-l 


where: 


F,  -  total  energy  of  the  frequency  band  at  the  source 
o 

k  -  coefficient  of  energy  dissipation 

A  -  Fourier  amplitude 

A  -  epicentral  distance  in  radians 

-  constants  for  each  frequency 

Taking  natural  logarithms  of  equation  (A-2) 

In  Ec  -  Dk  -  In  C,  +  2  In  A  +  In  sin  A  (A- 3) 

2 

Letting  E  =  C,A  and  C_  2lnA  ,  then 
o  1  o  2  o 

2lnA  4-  InsinA  =  C^-Pk  (A-4) 


A  least  squares  fit  of  the  points  obtained  at  several  azimuths  and  distances 
will  yield  an  average  k  (attenuation  coefficient)  for  the  region. 


TWO- STATION  METHOP 


Py  simply  equalizing  the  amplitude  between  two  stations  (Tsai 
and  Aki,  1970a)  we  obtain 


0(0 


7rf(r2  "  rl' 
T'(f) 


/ 


In 


A^f)  ^sinA^ 

TTTfT  V  sin^2 


(A-5) 


where: 


U(f) 

At(f) 

A2m 

rl’  r2 
Af  A2 


group  velocity 

spectral  amplitude  at  station  l 
spectral  amplitude  at  station  2 
epicentral  distances  in  km 
epicentral  distances  in  radians 


This  method  demands  high  quality  spectra.  It  is  also  highly 
desirable  to  have  several  stations  along  the  great  circle  path  to  eliminate 
minor  spectral  perturbations. 
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APPENDIX  B 


DETERMINATION  OF  FUNDAMENTAL  MODE  AMPLITUDE  SPECTRA 
AND  GROUP  VELOCITIES  BY  NARROW  BAND  FILTERING 


It  has  been  shown  by  Evernden  (1953)  that  the  propagation  paths 
for  surface  waves  are  greatly  altered  by  lateral  velocity  changes  along  the  path 
and  depend  on  frequency.  Further  evidence  of  such  effects  was  given  by 
Alexander  (1963),  Glover  and  Alexander  (1969),  and  Newton  (1973).  The  result 
of  this  phenomenon,  known  as  multipathing ,  should  be  analyzed  as  a  function 
of  frequency.  The  effect  of  multipathing  is  usually  an  introduction  of  anomalies 
(short  period  peaks)  into  the  power  spectrum  of  the  signal.  Because  the 
determination  of  source  parameters  from  far-field  radiation  patterns  involves 
studying  fundamental  mode  amplitude  spectra,  a  method  for  de- multipathing 

a  spectrum  is  required. 

One  method  which  has  proved  to  be  successful  in  de-multipathing 
is  the  application  of  a  series  of  narrow  band  filters  [Alexander  (1963), 
Archambeau  et  al  (1965),  Newton  (1973)].  Given  a  signal's  time  series  f(t) 
and  power  spectrum  F(oj),  one  chooses  a  frequency  to  be  examined  wq.  A 

narrow  band, 


“-[-a--4?-  -o+^] 

defines  a  unit  rectangular  function  by 

H(a>)  =  1  if  ojC  B  [then  F(w)  =  F(co)  H(oj)] 
=  0  otherw  ise. 


B-l 


If  Aoj  is  sufficiently  small,  the  inverse  Fourier  transform  of  F(oo)  will  be  a 

A 

well  behaved  oscillating  function,  f (t) ,  which  attains  one  or  more  distinct 
local  maxima 


fl  =  f(V-f2  =  f(V  -  f3  =  ‘  *  * 

N 

The  times  [tj.  indicate  the  relative  arrival  time  of  wave  groups  represented 

by  the  superposition  of  frequencies  in  the  band  B.  When  multipathing  occurs, 

2  7 T 

N>1  for  some  periods(s),  Tq  s  • 

o 

Now  suppose  one  chooses  as  the  frequencies  to  be  examined 

w  =  T**1  .  Zir  =  (10  +  5K)"1  .  2  *  K  =  l,  2,  ...  10 

o,K  o,K 


(or  periods  (T)  10  to  60  seconds),  and  we  want  the  wave  packet  functions 

f  ,  to  exhibit  M  cycles.  From  Fourier  transform  theory, 
l,  K 


At  Aw  «  47 r 
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A  oo  = 


At 


and 


o» 


M  =■ 


o,  K 


2  it 


At  = 


At 


o,  K 


A‘  =  To.KM 


Therefore 


A co  =  4 it  /T  K  =  l,  2,  .  .  .  10 

K  o,  K 


N„ 

Using  this  technique,  we  obtain  the  relative  arrival  times  [t  ]  ^  and  using 
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the  station  to  source  distance  (A),  the  origin  time  of  the  event  (tQ ) ,  the  trac 


start  time,  t  ,  we  obtain  group  velocities 
s 


N(K) 


T  a  ~rvi 

I  i,  K  “  t  -t  +t  .  . 
L  soi,  KJi=l 


K=l,  2,  ...  10 


Comparing  these  group  velocity  points  to  average  group 

velocity  curves  enables  one  to  choose  a  "best"  U.  and  corresponding 
7  i.  K 

A 

f.  „  for  each  frequency  (K). 
i,  K 


APPENDIX  C 


A  DISLOCATION  EARTHQUAKE  SOURCE-TIME  FUNCTION 


If  we  designate  the  source-time  function  of  the  theoretical 
earthquake  dislocation  as  D(t),  then  the  time  function  of  the  theoretical  force 
system  (single,  couple,  douple  couple,  etc.  forces)  which  causes  that 
dislocation  should  be  proportional  to  the  second  derivative  of  D(t);  therefore 


F(t)or 


d2D(t) 

dt2 


(C-l) 


Thus,  if  a  unit-step  function  (which  is  the  limiting  form  of  a  ramp  function) 
is  chosen  for  the  theoretical  dislocation,  i.  e. 


D(t)  =  CjlRt)  *  Jim  C.  1-e 

T-*  0  l' 

=  0 


t  >  0 
t  <  0 


(C-2 ) 


where:  C^  and  r  are  constants 

U(t)  =  1  t  >  0 

=  0  t  <  0 


then  the  time  function  of  the  force  system  will  be  a  doublet  function  according 
to  (C-l).  Therefore, 

F(ti « 

at2 


a 


dfi(t) 

dt 


ftp  (t) 


(C-3) 


C-l 


where:  8(t)  is  a  Dirac  delta  function 


(t)  is  a  doublet  function 

Now,  let  the  spectral  amplitudes  of  the  surface  waves  due  to 

the  douple  couple  force  of  1  dyne-cm  seismic  moment  with  temporal  variation 

of  a  Dirac  delta  function  as  given  by  Saito  (1967)  or  Harkrider  (1964)  be 
DC 

Ag  (w).  Then  the  spectral  amplitude  of  the  surface  waves  due  to  double  couple 
forces  with  general  temporal  variation  as  Ly(t)  will  be 


(to)  =  My  Ly  (to)AgC(w) 


Ijcn 


where  Ly (to)  = 
with  Ly  (t)  =  My  L y  (t) 


(t)e  dt 


(C-4) 


For  douple  couple  forces  varying  with  time  as  a  doublet 
function  as  given  in  (C-3),  corresponding  to  a  unit-step  function  of  the 
dislocation,  then  according  to  (C-4),  we  have 


A  (w)  =  (M*  w  )  Ac  (oj) 

My  /  ° 


(C-5) 


since  Ly(t)  -  Myft^(t),  then  Ly  (to)  =  oj, 
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APPENDIX  D 


REPRESENTATIVE  EARTH  STRUCTURES 

In  order  to  determine  the  effect  of  earth  structure  on  the 
source  parameters  obtained  from  either  an  exhaustive  or  non-linear  regression 
approach,  a  series  of  representative  structures  must  be  analyzed.  Because 
we  are  mainly  concerned  with  events  originating  in  Eurasia,  four  continental 
structures  have  been  chosen. 

A  Gutenburg- Bullen  standard  earth  model  is  shown  in  Figure 
D-l,  with  its  normal  high  velocity  crust.  A  thicker,  higher  velocity  crust 
on  top  of  this  standard  model  is  shown  in  Figure  D-2.  This  model  is  a  good 
approximation  for  the  Tibetian  Plateau.  Figure  D-3  shows  a  model  for  the 
Central  United  States  (McEvilly,  1963)  which  differs  from  the  standard  model  by 
higher  crustal  density  and  shear  wave  velocity  (V8).  Finally,  a  very  low 
velocity  crust  reflecting  sedimentary  layers  is  given  by  the  Hamiiton-Healy 
Earth  Model  in  Figure  D-4. 

Rayleigh  and  Love  group  velocity  curves  are  shown  in  Figures 
D-5  and  D-6,  respectively.  We  see  that  while  the  Central  U.  S.  and  Gutenburg- 
Bullen  structures  yield  almost  identical  curves,  the  Hamilton-Healy  yields 
lower  velocities  at  shorter  periods  while  the  Tibetian  Plateau  yields  lower 
velocities  at  longer  periods. 
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FIGURE  D-5 

THEORETICAL  RAYLEIGH  WAVE  GROUP  VELOCITY  CURVES 
FOR  VARIOUS  EARTH  STRUCTURES 
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FIGURE  D-6 

THEORETICAL  LOVE  WAVE  GROUP  VELOCITY  CURVES 
FOR  VARIOUS  EARTH  STRUCTURES 


appendix  e 

REVIEW  of  TSAI'S  METHOD 


Tsai's  method  for  the  estimation  of  seismic  source  parameters 
namely  focal  depth  00.  dip  angle  a.ip  angle  strike  and  seismic 

moment  (Ml,  is  the  exhausted  searching  scheme  in  finding  the  best  fit 

;rr  °bSerVed  3nd  "  ^  wave  amplitude 

pectra  for  one  or  more  observation  sites  TU  ,  .  , 

sense  that  .  °n  sites.  The  exhausted  search  is  in  the 

ense  that  every  combination  of  £  \  a  ,  ,  . 

nation  of  fi,  A,  cf>,  and  h  in  a  predefined  set  of  these 

;;;;;:ters  has  b~n  -  - «-  *  -  of 


et  of 

8, 

</>,  and  h  is 

usually 

defined 

by 

V 

60°  -  90° 

~ 

10° 

k  = 

4 

- 

90°  -  .90° 

= 

w 

o 

o 

!  = 

7 

V 

0 

-  180° 

A(h  ~- 
n 

10° 

n  = 

19 

3d 

3 

II 

0 

-125  .’an 

Ah  = 
m 

5  km 

m  = 

2i 

hus,  there  are  13832  members  in  the  set  A  ,wa-i  j  - 

he  S6t-  A  det*iled  description  of  Tsai's 
method  is  given  by  Tsai  and  Shen  (1972  b>. 

,  "y  TSai’S  me,h0d’  B“  -timation  of  seismic  source 

parameters  will  be  obtained  by  examining  28  e  which  have  h  „  • 

through  the  least-squares  fitting  „roc  • 

8  prncess’  Picking  up  the  Minfe,  ,1  and 

assigning  5k,  y  y  hm  and  M  the  values  which  correspond  to  that 

'k/!'  Th“  =  .  it  is  strictly  a  minimum  s  criterion.  The  distribution 
and  magnitude  of  these  28  «fc/  have  no.  been  taken  into  consideration 
ecause  the  distribution  is  incomplete  and  the  magnitudes  of  the  28  <  are 
Widely  scattered  with  little  trend  to  follow.  kf 
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k  t 


In  order  to  attain  a  greater  degree  of  confidence  in  using  the 

minimum  residual  Min  ft  ^  solution,  a  detailed  knowledge  of  the  residual 

dist  ribution  would  be  helpful  in  determining  the  distribution  about  the 

absolute  minimum  and  guard  against  a  possibie  spurious  minimum  generated 

by  poor  quality  data,  Examining  only  those  28  t  which  correspond  to  the 

predetermined  28  combinations  of  dip  and  slip  angles  somehow  impHes 

beforehand  knowledge  obtained  from  a  body-wave  fault  plane  solution.  ,t 

not  affect  the  optimal  source  parameter  estimation  based  on  the  Min(f 

enterion.  but  it  will  definitely  make  differences  in  the  distribution  and 

magnitude  of  .fc/  .  Tor  events  originating  in  Eurasia,  body-wave  data  with 

good  azimuthal  coverage  are  usually  no,  immediately  available.  Therefore 

,f  W°  are  ',cee"di"»  °"'V  -  -Taco  wave  data,  residua,  distribution  informa 

,l°"  ‘S  .imP°r,anl-  THe  changes  were  instituted  in  Tsai's  method 

to  obtain  optimal  source  parameter  estimation: 

•  In  the  fitting  process,  the  focal  depth  will  be  varied  last, 

because  the  depth  <h)  will  effect  the  shape  of  the  theoretical 
surface  wave  amplitude  spectra  the  most,  especially  for 
Rayleigh  Waves  (see  Section  III).  There  will  be  one  Ch 
for  each  group  at  one  focal  depth  (there  will  be  26  groi^' 
with  532  members  in  each  group  for  the  set  of  fi,  A,  <f>,  and 
h  described  earlier). 

Ihe  variation  of  with  respect  to  focal  depth  will  be 

obtained  and  a  standard  value  of  e  reh  .  1  ... 

,  lf  h  ’  '  min  'standard,  wlU 

bo  determined  from  those  e  fnr  . 

muse  e  min  '°r  the  purpose  of  picking 
appropriate  members  in  the  set  of  5,  A,  0,  and  h  to  form  a 
sample  group, 

•  Members  in  the  set  of  6,  A.  d>.  and  h  whose  (  is  less  than 

(  min'sTANnARO  wl11  be  t,ktn  int<>  consideration  in  forming 

the  distribution  of  t  with  respect  to  focal  depth. 

Examples  of  these  distributions  for  each  source  parameter  are 
Riven  in  Sections  IV,  V,  and  VI. 
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APPENDIX  F 


NON-IJNEAR  REGRESSION  OF  SEISMIC  SOURCE  PARAMETERS 
FROM  SURFACE  WAVE  DATA 


We  will  assume  that  the  seismic  source  as  observed  in  the 
far-field  can  be  represented  by  a  simple  double  couple  in  the  form  developed 
by  Harkrider  (1964).  This  representation  may  not  be  the  most  accurate, 
but  for  this  discussion  we  will  assume  that  higher  order  components  are 
not  needed  (Turnbull,  1972).  The  general  form  for  the  far-field  Rayleigh 
and  Love  wave  radiation  from  a  double  couple  source  may  be  written  as: 

.  3rr  -iknr 

U  =  Skne"lT  N(h)X(fl,h)- - —  (F-I) 

r  2 


whe  re 

U  :  spectral  horizontal  Rayleigh  or  Love  wave  displacement 
S  :  spectral  source  function  (For  step  function  -  M /io ) 
k^:  (,J/c  Rayleigh  or  Love  wave  angular  wave  number 
N(h)  :  either  N0(h)  or  N  r  z  (li )  -  response  functions  from  the  layered 
half  space  for  Love  and  Rayleigh  waves  respectively 
X(0,h)  :  the  complex  radiation  pattern  function  which  can  be  written  as 


X(0,h)--d  +  i(d.sin0  +  d  cos  0)  +  d.sin  2  0  +  d  .cos  2  0  (F-2) 

o  1  L  3  4 
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The  coefficients  d.  are  given  below  for  a  double  couple: 


d 

o 

LOVE 

RAYLEIGH 

0 

sin  A  sin  28B(h) 

d! 

cos  Acos8G(h) 

-  sinXcos  28C  (h) 

d? 

-  sinXcos  28G(h) 

-cosAcos8C(h) 

d. 

sinXsin  2  8V (h) 

cosAsin8A(h) 

d4 

cosXsin8V(h) 

-■£  sinXsin  28A(h) 

where  d  is  the  azimuth  angle,  8  the  dip  angle,  X  the  slip  angle,  and 


A(h)  = 

R(h)  = 

C(h)  = 

V(h)  = 

G(h)  = 


-[\  <h)/wo] 

(3-4|)^,wH+iH8(h,/WCr5 

[^sm)/(w0/crj 

[v8(h)/v0] 

(’ 


|vo/Cl 


where: 


u  (h)  w  (h)  v  (h) 

8  S  S 


w 


w 


-  Haskell's  plane  wave  particle 
velocity  (or  displacement)  ratios 


ar  Ih)  -  Depth- dependent  factor  of  normal  stress  associated  with 
the  Ravleieh  waves 

7  (h),  r,  (h)  -  Depth-dependent  factors  of  tangential  stresses 

i  c  Xq 

associated  with  the  Rayleigh  and  Love  waves. 
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The  fault  plane  geometry  is  shown  below. 


FAULT  PLANE  GEOMETRY 
After  Harkrider  (1964) 

If  we  are  to  fit  the  radiation  patterns  described  by  the  displace¬ 
ment  function  U,  we  must  determine  which  parameters  would  be  feasible  to 
use  in  a  regression  analysis.  We  first  realize  that  a  fit  must  be  made  for 
each  frequency  (or  wave  number),  because  radiation  patterns  vary  with 
frequency  for  an  earthquake.  [This  is  due  to  the  effect  of  rupture  velocity 
(Archambeau,  1964)].  Also,  the  variation  of  depth  in  the  regression  procedure 
requires  a  significant  effort,  because  a  different  placement  of  the  source  in 
the  layered  half  space  would  involve  a  complicated  variation  of  the  response 
functions  N(h),  A(h),  R(h),  C(h),  G(h),  and  V(h).  As  a  first  attempt,  the 
depth  is  estimated  using  an  exhaustive  technique  such  as  Tsai's  method. 
Regression  is  then  performed  at  each  of  several  reasonably  spaced  depths 
bracketing  this  estimated  depth,  the  final  depth  being  that  which  produces  the 
minimum  residual. 
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The  parameters  which  easily  lend  themselves  to  this  type  of 
repression  analysis  are  A(slip  angle),  §(dir'  angle),  ^(strike  angle),  and 
V(seismic  moment).  TTpon  examination  of  expression  (F-2),  we  see  that  a 
linear  regression  approach  would  involve  significant  errors  because  the  angles 
are  arguments  of  trigonometric  functions. 

There  are  several  approaches  available  for  non-linear  regression, 
but  possibly  the  most  direct  and  easily  understood  involves  making  a  Taylor 
Series  approximation  ^Draper  and  Smith,  1968).  The  objective  of  this  procedure 
is  to  obtain  the  expressions  involving  strike,  dip,  slip,  and  seismic  moment 
using  the  first  terms  in  Taylor  Series  expansion,  and  obtain  values  for  the 
regression  parameters  with  an  iterative  technique. 

Py  assuming  the  theoretical  description  of  the  source,  we  have 
the  set  of  theoretical  equations 

v  =  f'x.,  /3j)  (F-3) 


wh  ere: 


x  =  x,,  x  ....  x  independent  variables  (in  our  case,  we  have  depth 
i  1  2  n 

and  azimuth). 


/?j  Py  P2> 

and  slip). 


.  p  regression  parameters  (in  our  case,  strike,  dip, 


The  objective  is  to  find  the  estimates  of  which  minimize  the 


n  n  ^ 

error  of  squares"  =  S(/T)  =  ^  Fu  =  ^  (Yu  -  fu^Xi»  P ;))  (F-4> 

u=l  u=l 


where: 
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-  observed  dependent  variables  (u=l,  2,  . 


•  9 


n)  (in  our  case,  the 


v 


amplitude  variation  with  azimuth  for  a  particular  period) 


F  -  observed  deviation  about  predicted  equation 
u 

vu  -  /3j)  +  for  each  u  (u  -  1,  2,  .  .  .  ,  n) 


(  v-  5 ) 


Regression  satisfying  condition  (F-4)  may  be  achieved  using  the  'linearization' 
or  'Taylor  Series'  method. 


Carrying  out  a  Taylor  Series  expansion  of  f(x.,/3.)  about  the 

0  J 
point  and  curtailing  expansion  at  the  first  derivative  yields 


n  £ 

'u'Vft'  ft'Vft'  +  L, 

j=l 


dP\ 


<ft-  ft°> 

o 


ft  ■  ft 


0 

where  ft.  and  /J.  are  separated  by  the  small  increment  A. 
Now  setting 


f°  -  f  (x.,/3°) 
u  u  1  'j 


(F-6) 


(F-7' 


A.  -  B.  -  B?  (  F-8) 

J  J  'J 


0 

7. 

J  • u 


ft 'ft 


(F-9) 


yields 


Y 

u 


u 


(  F-10) 
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which  is  a  'linearized  form  of  expression  (F-5). 

Applying  the  least-squares  condition  (F-4)  to  (F-10)  gives  as  the  "sum  of 


squares 


"  r  o  ^  o 
ss</y ’E  vv£vzi, 

J  u-l  J=1 


(  r  -H) 


Minimizing  F  with  respect  to  A  gives 


i.2f  fY  r° _  y  a  •  z°  1  [y;  z° 

m  u  »  f_ i  j  >uJ  [a  j’u_ 


(F-12) 


£Mi-£  i‘,<. 

u- 1  L  *'  u=l  j  =  l 


(F-13) 


as  a  no 


n-trivial  solution  of  (F-12).  These  expressions  (F-13)  are  sometimes 


called  the  normal  equations. 

For  purposes  of  clarity  and  commutation,  we  can  write  equation 
(F-13i  in  matrix  notation.  Therefore,  we  have 


01  " 

Pz  -P°2 

a  P° 

pP  - 

'  (pxl)  ' 


(F  - 14) 


p  (regression 
parameters) 


F-6 


In  our  case,  p  =  4  for  <f>.  \,  8,  M. 


-  P$ 
h  -  fix 

h  "  P& 

f)  '  ° 

‘  Pm 
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0 


(regression  parameters) 
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0 

11 
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z° 


in 


Z 


...  7° 

21  pi 

. 7° 

12  p2 

. 7° 

pn 

(nxp)  or  (nx4)  — 1 


n  (number  of 
observations) 


( F —15) 


In  our  case, 
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0fn 


Zn 


ff[n 


80 


A 


~  A»  8.  M 


P.  ~  <t>,  6.  M 
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3n 
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4n 


£fn 

afn 


Pj  =  <t>,  A,  M 


/J.  ~  A,  8.0 
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Y  - 


‘l 

f2° 

.0 


i—  n  n  — 1 

(nxl) 


t 

n  (number  of 
observations) 


(P  -16) 


F-7 


Then  expression  [F-131  is  given  by 


(F-17) 


The  above  equation  may  be  solved  iteratively  for  estimates 

,0 


/?•  ot  the  true  regression  parameters  p.  by  making  a  first  guess  at  p  .  a 

.  J  o  J 

using  (F-17)  to  give  fl..  This  is  used  as  a  new  estimator  R.  which  is  us 

J  H J 

calculate  p.  and  so  on  until  the  iteration  converges. 


nd 

ed  to 


Before  we  proceed  further  with  this  discussion,  a  note  must  be 
made  on  applying  this  technique  to  our  source  description.  We  see  from 
equation  (F-2'  that  the  radiation  pattern  is  given  by  a  complex  expression. 

This  form  does  not  lend  itself  to  a  regression  analysis.  Instead,  we  can  work 
with  the  modulus  squared.  Therefore,  we  are  working  with  a  quantity  directly 
proportional  to  the  energy  in  the  seismic  wave  rather  than  the  amplitude. 

We  can  therefore  write 


X(fl)X  (ft)  -  d  +  d  sin  26  f  d  cos  2 ft  4-2dd  sin20  +  2d  d  cos20 
”  j  4  3  0  4  0 

*■  2d3(*4  sin  2 0  cos  2ft  4  dj  sin  2ft  +  cos  26  +  2djd2  sin0cos0 


( F- 18) 


with  the  d.  defined  as  before.  Then  X(ft)  X  ‘  (0)  will  be  the  function  f  (x.  ,R  ). 

L  u  x  C’i 
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APPENDIX  G 


REPRESENTATIVE  SPECTRA  AND  RADIATION  PATTERNS 

The  following  spectra  and  radiation  patterns  were  generated 
using  a  theoretical  double  couple  point  source  placed  in  a  Guttenburg  layered 
half  space: 

Figures  G-la,  b,  c  to  C-4a,b,  c  -  Rayleigh,  Love,  and  Love/Rayleigh 

amplitude  spectra  for  a  vertical  strike- slip  fault  at  30  azimuth 
showing  variation  of  depth,  dip,  slip,  and  strike  respectively. 

Figures  C-5a,  b,  c-  Rayleigh,  Love,  and  Love/Rayleigh  amplitude 
spectra  for  an  approximately  vertical  strike-slip  fault  at 
30°  azimuth  showing  variation  of  depth. 

Figures  G-f.a,  b,  c  to  G-9a,  b,  c  -  Rayleigh,  Love,  and  Love/Rayleigh 
amplitude  spectra  for  a  dip-slip  type  fault  at  30  azimuth 
showing  variation  of  depth,  dip.  slip,  and  strike  respectively. 

Figures  G-lOa,  b  to  G-I4a,  b  -  Rayleigh  and  Love  wave  radiation 

patterns  for  h  =  15  km,  ~  f>0  ,  A  -  -30  ,  and  f/>-  20  for 
periods  50,  40,  30,  20,  and  10  seconds  respectivelv. 
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FIGURE  G-lb 

LOVE  WAVE  AMPLITUDE  SPECTRA 
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LOVE/ RAYLEIGH  WAVE  AMPLITUDE  SPECTRA 
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FIGURE  G-5c 

LOVE/  R  AYI.EIGH  WAVE  AMPLITUDE  SPECTRA 
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FIGURE  G-6b 

LOVE  WAVE  AMPLITUDE  SPECTRA 
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FIGURE  G-8a 

RAYLEIGH  WAVE  AMPLITUDE  SPECTRA 
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FIGURE  G-8c 

LOVE/ RAYLEIGH  WAVE  AMPLITUDE  SPECTRA 
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FIGURE  G-9a 

RAYLEIGH  WAVF  AMPLITUDE  SPECTRA 
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FIGURE  G-IOa 

RAYLEIGH  WAVE  RADIATION  PATTERNS 
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FIGURE  G-tla 

RAYLEIGH  WAVE  RADIATION  PATTERNS 
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FIGURE  G-Ub 


LOVE  WAVE  RADIATION  PATTERNS 
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FIGURE  G-  12a 

RAYLEIGH  WAVE  RADIATION  PATTERNS 
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FIGURE  G-13a 
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FIGURE  G- 1 3b 

LOVE  WAVE  RADIATION  PATTERNS 
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APPENDIX  H 


EXAMPLES  OF  MU LT IP ATHING  FROM  SINKIANG 
TO  SELECTED  VLPE  STATIONS 

For  practically  all  VLPE  records  of  Sinkiang  events,  severe 
multipathing  occurred.  Figures  H-l  show  the  results  of  narrow  band  filters 
applied  to  the  record  at  CHG  from  event  KIRSI/059,  with  Figure  H-2  showing 
the  group  velocities  of  the  wave  packet  peaks  at  each  period  plotted  with  stan¬ 
dard  group  velocity  curves.  The  wave  packet  peaks  are  numbered  in  order  of 
descending  amplitude.  We  see  that  by  matching  a  general  trend  of  the  curves 
with  the  plotted  points,  direct  path  amplitudes  are  obtained.  For  both  the 
Rayleigh  and  Love  waves  (Figures  H-3:  narrowband  filters,  and  Figure  H-4: 
Love  wave  group  velocities)  recorded  at  CHG,  we  see  very  low  group  veloci¬ 
ties  (although  some  of  this  could  be  due  to  group  delay  of  the  instrument).  Two 
other  sets  of  Rayleigh  and  Love  wave  group  velocity  curves  are  shown  in  Fig¬ 
ures  H-5  to  H-8,  demonstrating  'typ  paths  encountered.  Figures  H-5  and 
H-6  for  event  KIRSI/059  recorded  at  FBK,  show  fairly  clos  i  adherence  to  the 
Rayleigh  wave  group  velocity  curve,  while  severe  multipathing  occurs  for  the 
Love  waves  below  30  seconds  period.  Figures  H-7  and  H-8  for  the  same  event 
recorded  at  KIP  show  the  effects  that  a  part  continental-part  oceanic  path  have 
on  the  group  velocity  curve. 
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FIGURE  H-lb 

NARROWBAND  FILTER  ANALYSIS  OF  THE  VERTICAL  RAYLEIGH  WAVE 
OF  EVENT  LX/KIRSI/059  RECORDED  AT  CHG 
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